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PRESSURE DISTRIBUTION ON JCUKOWSKI WINGS.* 
By Otto Blumenthal. 

In the winter semester of 1911-12, I described, in a lec- 
ture on tue hydrodynamic bases of the problem of flight, the 
potential flow about a Joukowski wing.** In connection with, 
this lecture, Karl Toepier and Erich Trefftz computed the 
pressure distribution on several typical wings and plotted 
their results. I now publish these diagrams accompanied by a 
qualitative discussion of the pressure distribution, which 
sufficiently indicates the various possible phenomena. For a 
quicker survey, I have divided the article into two parts, 
the first part dealing with the more mathematical and hydrod^r- 
namic aspects and the second part, which is comprehensible in 
itself, taking up the real discussion from the practical sta.nd- 
point-. 



* From "Zoitschrift fur Flugtechnik und Motorluf tschif fahrt. 
May 31, 1913. 
** See above magazine. Vol. I (1910), p. 281- 



P.A^Ci.A. Tochnicc^l McinoraiKium No* 336 



2 



I 

V/e" obtain the entire number of all Joulcowski wings of the 
length 3 1 with the trailing edge at the point x = - I, by 
laying^ in a ^ = E + i T} plane through the point ^ = l/2, 
the cluster of all the circles vvhieh contain the point X = 
either inside or on their circumference, and plotting these 
circles. by nieans of the formula 

on the z = X + i y plane. The circles, which contain the 
point ^ = 1/2 on their circtomf er ence, thus become doubly in- 
tersected 8.rcs and, in particular, the circle, which has the 
distance (-1/2, + l/Z) for its diameter, becomes the recti- 
linear distance of the length 21 . The circles which contain 
the point ^ = l/2 inside, furnish the real Joukowski fig- 
ures. The point ^ = - l/2 passes every time into the sharp 
trailing edge* The individual Joukowski wings are character- 
ized by the following quantities (Fig. l). The center M of 
the circle K is connected with the point H, ^ - - l/2, and 
the point of intersection ,of this connecting line vrith the rj 
axis is designated by . The distance OM' on the il axis 
is equal to half the height of the arc produced by describing 
the circle about M' as its center and is therefore designated 
by f/2, as half the camber of the Joukowski wing, f being 
its first cha>ract eristic dimension. We have chosen as the 
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second characteristic dimension, the radii difference im^ = 6. 
This gives a measurement for the thickness of the Joukowski 



wing. 



7/e will now consider the determination of the velocity 
and pressure distribution which produce an air flow along the 
wing, in infinity, with the velocity V at an angle of tr-p 
with the positive x axis, g being the angle of attack of the 
wing. 

The absolute velocity q of this flow is calculated thus: 
If ^- {i ,'']) is the absolute velocity of the air flow, of ve- 
locity V and angle of attack around the circle K in 
the C plane, then 

q (x,y) = ^ ' • 



az 



It is, however, 
ds 



and a.ldng the circle K 



K a ,r,) = 



2 V (I sinp + T] cosP) + c|. 



+ 6 



where 3 rr c is the circulation. This constant is dotcrrained 
according to Kutta, by the condition that the- velocity at the 
trailing edge is finite and therefore, since dz/dC there dis- 
appes-rs, k must also disappear at the point H. Thus we ob- 
tain 
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K (e,r,) = 



- . i ZU + ^) sinP + Ti cosp 



q 

V 



a' 



& + ^ ; sinp + ri cospl 

i 



\ 4 / 



2 



g2 . ^3 



(3) 



This rowthor involvccL expression is simplified by the in- 
troduction of a new variable, the angle 00 at the center of 
the circle measured from the radius ISR. In this angle, 

the coordinates i and the quantities connected with them 
are expressed as follov/s: For abbreviation, v/e designate the 
radius of the circle K \Yith r == — + -S and introduce 



the angle A by 



cos A = 



sin A 



The geometric significance of A and o) is obvious from Fig. 1- 
By simple calculations we now obtain 



l=-|.2rsin|sln(|.A), 



T) = - 2 r sin g cos + a\ 



(3) 



0^ = 21 r sin 3ini''^- + + 4 sin^ ^ 

4 2 \(': / 2 



— + 4 r sm - [ S sm - - - cos ( - + A J 



(4) 
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Formula 3 for q is simplified by the introduction of the 
angle to 



3 











V 2 / 


— — — ' ' — "cr- 



COS' 



6 sm 2 - 2 



4 



3 * a) 



(5) 



From this v;e next derive a few general results which hold 
good for all the quantities f , 6 , p . 

a) On the trailing edge £ = ^ cos (A + P). 

V ^r 

b) On top of the wing, there is always a portion along 
which the velocity q > V, hence where there is a negative 
pressure* As proof of this, we will consider the center of 
the upper side, the poixit o) = ^ - A. At this point a > r. 



a 



s can easily be seen geometrically (Fig. l) or from formula 4. 



However, if we put ^ ~- - k in formula 5, it then becomes 



q 

V 



C0£ 



r i 



oc; CH. - k\ At . '4 ^ f>f 

°^ V A 2 / y 1 + -55- ( 0 + - ) 



2/ 



COS - ^ - 6^ 



COS ' 



/TT 



\ 4 



3 / 



COS 



A 
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TT 



COS 
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c) The velocity is zero at the point o) = n - 2A - 
which is always located on the lower side. At this point the. 
streamline enters the v/ing. Further general conclusions (i.e., 
applying to all f, 6, |3 ) can hardly be drawn* V/e obtain 
considerably more accurate expressions in the especially int- 
eresting practical case where^ in the vicinity of the leading 
edge J a pronounced velocity maximum and consequently a strong 
suction is produced. We vrill confine ourselves to this case 
in all tlia,t follows. Hereby we can, in formula 5, first of 
all disregard the slight fluctuation of the factor for 
small values of f and. 6 and consider only the factor 
vfhich imist be alone decisive for the great changes in velocity?-. 
This factor, however, enables a simple explanation. 

For this purpose, we introduce the angle == |- + A + P . 
The entering point of the streamline then lies at t = '^/S, 
v/here F disappears- In general, we have 

\- = (a^ + sin^fB) tan^ t - 2 (ab - sinp cosP) tan^k + 
F 

+ (b^ + cos^ P ), 

( 6^ 

a = ^ f 6 COB {A B) - ^ sin{3 \, ^ . ) 

> 

b = - 6 sin (A -i- B) + t cosp^ . 

Consequently, F attains its m.aximum value at the angle 
which is given by the formula 

tan % = ab^^inB cosp (7) 
a^ + s in^ j3 
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and this value is 



/ (5 cos (A + ") - I sinB^ + -t- sin^p 

F-:r.ax = , - ' ' ( 7' ) 

. 5 cos A 



We now make the assunption, corresponding to the already 
announced purpose of our investigation, that F has a hig;h 
mximtim in relation to the value of cos (A + on the trail- 
ing edge* We require, e.g., that Fp^^x shall equal or exceed 
y~2. This is inathematioaj.ly the most favorable. Formula 7', 
with the c-id of e. rough estima.te, then gives 



Bine ^ 6 cos (A + 3) fl + ——1=^] 

= 5 cos (A + P ) (14-- sin A) .... (8) 
V/ith this insertion, the nurnoratpr. of tan ij/^ is smaller tha.n 

^ 46 



^ cos (A + P ) 



•LJ^^±^1 cosP - 6 sin (a + P ) (l-sinA) 



For small f, 5 and B, this value is always negative and there- 
fore the maximm value of F is assumed to be at a point lo- 
cated between the entering point and the trailing edge on the 

portion of the sur face belonging to the upper side.* On the 

* Generally the point is located on the upper side. It lies 

between the entering point and the leading edge, only when ' 6 
is very sma.ll in comparison with f. For 6 = 0, it lies on 
the leading edge. 



NsA.C.A. Technical Hcrnorandum No. 356 



8 



other hand^ it can be shown that the maximum is located not 
far from the entering point- In fact the greatly pr epond.crat- 
ing mcmiDer in the numerator of tan \^/^, on account of formula 
8, is cos3 sinfi • The case is not quite so simple with the 
denominator, which is 

a^ + sinPp 

= ^ S'' cos^ (A + - 6 f cos (A + p)sin!3 + ^ ^ sin^p . 
U^* L 4 J 

If we introduce into the first mem.her, on the right side of 
formula 8, the above limit for 6, the denominator is then 
smaller than 2 ^ sin^ P • Hence tan is either siraller . 

or at most only unessentially* grea^ter than - ^ cotp, . which 
shows that t is cither smallc.r or at most- only slightly great- 
er tkan ^ + 2 p. The point (o , at which F assumes its max- 
imum value, is located between the entering point of the stream- 
line and the upper side and, at most^ only slightly farther 
than 4p from the entering point. 

Lastly, it m.ay be rema^rked that in formula 6 for r;^^ both 
the powers, tan^^ and tanx]/, appear to be multiplied by 
small coefficients. Both these members therefore assume quite 
large values for large values of tani, i.e., in the immedi- 
ate vicinity of the entering point and the leading edge. Hence 
F differs but little over the whole surface, with the excep- 
tion of the specified region, from the value cos (A + on 

the trailing edge and th erefore only slightly fro m unity . 

* "Unessentially" means that the deviations are of the order of 
magnitude f/l , 5/l . 



i 
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The results of the conputatione in I arc as follows: A 
Joukowski v;ing is characterized by the three dinensions, name- 
ly, the length 21, the camber f and the radii difference 
6, v/hich is expressed in the thickness of the v/ing. To these 
is added the angle of attack p . The points on the surface 
are most conveniently coniputed with the aid of a variable cjo of 
the an^le at the center of the circle in Fig. 1. The fornulas 
for the coordinates will not be given here* In practice, a 
graphic process is employed which is explained in the accompa- 
nying note by E. Treffts (pages 130-131 of this same volume 
of "Zeitschrift fiir Flugtechnik und liotorluftschiffahrt • " ) 
We require only the follov/ing data: w = 0 gives the trailing 
edge; 0) = TT - 2A (tan A = y) gives the leading edge and the 
intermediate values of o) correspond to the lower surface of 
the viTingo ~ tt 2A - 2 fB gives the entering point, i.e., 
the point where the air flow strikes the surface and hence 
where the velocity is zero. 

The ratio of the absolute velocity q of the air. flow 
on the wing to the velocity Y in infinity is given by 



X = ^' ^ "^^^^ radius of the circle in Fig. 1 and 



q 

V 




F 



(5) 



2 
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a the distance OB in the same figure. Both factors, a and 

F, depend on o) ^ For small f and 6, the value of the factor' 

— j- differs but little from unity; The propertiee of the fac- 
r — 
2 

tor F, as obtained by the calculations of I, can be summarized 
as follov/s. 

On the trailing edge, F has the value cos (A + p) and 
decreases at the customary angles of attack (about 6°), from 
the leading edge to the entering point, where it becomes zero. 
For small f and 6, the decrease takes place very slowly 
throughout most of the lower side and first becomes rapid in 
the immediate vicinity of the entering point. 

From the entering point, , F increases rapidly and attains 
near the leading edge, a maximum of the order of magnitude 
pl/26. This is approximately also the maximum value of the 
velocity ratio q : V. For this maxiraiim value, the ratio of 
the angle of attack to the thickness of the wing is therefore 
decisive, the camber having, in the first order, no effect on 
it. In constructive wing shapes, where 6 and P are of the 
same- order of magnitude, the velocity at the leading edge is 
accordingly not very great. This result is important because 
it explains the effect of rounding the leading edge* The re- 
suit is still more striking when. we consider the radius of 
curvature p of the leading edge- It is, namely, with unes- 
sential omissions £ = 16 -J fl - 4 4^. The radius of curvature 
therefore diminishes rapidly with decreasing 6, the rounding 
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off of the leading edge "being very slight, and the maximum ve- 
locity remains within moderate bounds. The negative pressure 
on the leading edge, which, according to Bernouilli's equation, 
is proportional to q^/V^, is computed by the introduction of 
the radius of curvature in the first approximation, to 
4 I consider this simple formula v/orthy of attention. 

The course of F along the top of the wing can- finally 
be characterized as follows: At some distance from the leading 
edge, F changes but slov^ly. If, therefore, the maximum value 
of F is much greater than unity, it falls abruptly at first 
and then gradually approaches the value at the trailing edge. 

Only in the vicinity of the leading edge does the factor 
F give us sufficiently accurate information concerning the 
course of the velocity q. Everywhere else we need to know 
ttte course of a . This can be easily found geometrically from 
Fig- 1. On the leading edge a lias the value of .1/2. From 
the triangle ffliO, it follows that a assum.es its minimum value 
for the angle o:- , which ds given by 



^min ^ A 



sin A '^^ o^cos" A + r" sin A 



Hence 



^min 



/ 4 X ^ £^ 



For the angle 2 ^rain^ '^'^ again have ^ = |r sl^^L then a 
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increases further, up to the angle Wjjj^^ = w^^^^^ + rr. The value 

^niin increases v;ith f/6 • For |- = 0, ta^^^ = 0, hence 
the value of a is smallest on the trailing edge and greatest 
on the leading edge. Conversely, for 0, '«>inin = - A 

and is therefore situated in the middle of the lower surface. 
In general, with increasing f/6, the minimum value of a 
moves from the trailing edge to the middle of the lower sur- 
face; the point where a - ^ again from the trailing edge to 
the leading edge, while the maximum value of a moves simul- 
taneously from the leading edge to the middle of the upper 
surface. 

In order to get an idea of the course of the velocity, we 
must now estimate the mutual effect of the factors F and a^. 
I v;ill proceed with this discussion in close connection with 
the diagrams, which I must first explain. Their arrangement 
is the same as for the diagrams in Eiffel's "Resistance de 
I'air." Each figure has, at the bottom, an accurate outline 
of the wing section. Vertically above each point of the wing, 
there is plotted from a zero line, on the vertical the ratio 
qVv^, the upper curve corresponding to the upper side and 
the lower curve to the lower side of the wing section. The 
dashed line shows the unit distance from the zero line. The 
area enclosed by the q^/V^ curve gives, when multiplied by 

7 2 

^ V , the lift of a unit width of the wing. The chosen an- 
gle of attack is 6°0 = O.l), the air flow being horizontal. 
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I am dividing tlac discussion into several paragraphs. 

1. The strong suction on the greater portion of the upper 
surface is couirnon to all the figures. This is indeed the chief 
source of the lift, while the pressure on the lov/er surface 
contricutes only a small increment. This can be easily veri- 
fied from the general lav/s* In fact, as already stated, F 
diminishes very slowly along the under surface from the trail- 
ing edge almost to the entering point. Hence. F differs but 
little, on most of the lower surface, from the value 

cos (a + P)^ which it has on the trailing edge, and therefore 

only a little from unity. Since also the factor falls 

r — 

2 

only slightly below 1, up to the vicinity of the entering ■ 
point, oyv is certainly not much smaller than one and hence 
there is only a slight pressure. 

6 1 

2. Fig. 2 (f = 0 and y = -j^) indeed shows a suction 
effect along a portion of the. under side. Since the values 
of f are here less than unity, such a suction effect ca,n 
only be very small. Its appearance is due to relatively large 
values of (f and depends essentially on the ratio f : 6. 
Between the trailing edge and w = 2(si^^^ no suction can occur, 
because a is here smaller than l/2. Any suction effect can 
therefore be expected for only small values of f : 6, where 
'^rnin small. For 6=0 any suction effect is entirely inv- • 
possible, since 2a)jj.^j^ then corresponds to the leading edge- 
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The suction effect has also been experimentally determined by 
Eiffel on the wing en aile d'oiseau" which probably alone of 
all the surfaces tested by him can be compared \vith a Joukorslci 
wing section.* 

3. Even on the upper side, the course of the velocity is 
characteristically affected by the ratio f : 5. This is clear- 
ly shown by Figs. 2-4. f/6 is expressed on the upper surface 
in the position of the angle ^-^.^x ~ ^min for which o 

has its ipaxinTum value. The fact that a continues to increa-se 
from the leading edge as far as ^-^.^^^x oauses the inaximum value 
of q to move farther from the leading edge than the maxintum 
of F and the fall in velocity to be less ra.pid» We differ- 
entiate "slightly cambered" wings (f/6<2),. in which co^^x 
near the leading edge, and "highly cambered" wings (f/5>3), 
in which ^xc^ax l^^s nearer the middle of the upper surface. 
Fig. 2, I'vith f = 0, is a typical example of a slightly cam- 
bered wing. Here the maximum value of a is situated in the 

leading edge and the V3.1ues of F and a therefore decrease 

* Eiffel, "Resistance do I'air," 1911, Table XII; alBo p. 105 
and "Complement," p. 192 ("aile Nieuport"). The Eiffel fig- 
ures show that the suction effect increases on the under side 
with decreasing angle of attack. This is in agreement with 
our theori-r, for the factor F increases, as shown by formula 
6, at every point on the under side with decreasing An- 
other suction effect, v;hich Eiffel finds on the trailing edge 
of nearly a,ll wings, is doubtless due to the formation of 
vortices. 
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simultaneously, thus producing a very pronounced maximum veloc- 
ity, a-1 though the maximum velocity is not important in itself. 
On the other hand. Fig. 4 (i/l l/5, 6/l = l/20) shov/s, in 
spite of a twice as large maximum velocity, a remrkably slow 
velocity decrease toward the upper side* In fact, in these 
experiments, the maximum of a is cituated at about l/3 of the 
upper side and a more rapid velocity decrease accordingly 
first begins behind this point. Wc note also the small inter- 
mec5.iate maximum, on the uxjpcr side, which is caused by the in- 
crease of in spite of the simultaneous decrease in F. The 
mean between Fig. 2 and Fig. 4 is held by Fig. 3, with 
f/l = l/lO and 6/l = l/20. Here ooniax does not lie very far 
from the leading edge, about l/5 of the upper side, the in- 
crease in vanishes under the decrease in F and along the 
greater portion of the upper side we note a uniform falling 
off in velocity, due to the simultaneous decrease in the fac- 
tors F 3-nd a^. 

These relations were also found in Eiffel^s experiments 
with the wing "en aile d'oiseau." Even the intermediate maxi- 
mum of q on highly cambered v/ings is found on his figures.* 

Lastly, I wish to call attention to the fact that the 
measurements of Fig. 4 appear to me to be worthy of commenda- 
tion, on account of the very uniform stressing of the upper 
side. 

* On Eiffel's figures, it appears that, with decreasing p, 
the mximum velocity moves backward from the leading edge on 
the upper side. This also agrees with the theoretical conclu- 
sions. 
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4. Fig. 5 has a very slight rounding (6/l = l/50) at 
f/l = 1/5. Therefore 25'^ 2.5 and hence the very high ve- 
locity inaximam on the leading edge.. Vve liave already seen that 
the high maxinia must decrease very rapidly toward the upper 
side. The region of this steep decline corresponds to an angle 
of about the size |3 . During the drop, however, there is in 
the figure aMong space of almost constant velocity. This is 
explained, as in paragraph 3, by the fact that the iTiaximum 
value of a is located at about 1/2.5 of the upper side. Only 
behind this point is there again a rapid decline to the trail- 
ing edge. This behavior is generally characteristic for high- 
ly cambered wings of slight rounding and occ\irs also on Eif- 
fel's diagrams. 

As regards the production of the diagrams, it may be noted, 
in conclusion, that they were drawn according to the very con- 
venient m.ethod of E. Trefftz, as set forth in the accompanying 
note* V/herever it appeared necessary, the plotting was veri- 
fied by calculation. 
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GRAPHIC CONSTRUCTION OF JOUKOWSKI V/INGS.* 
By S. Trefftz. 

In plotting the cross-sectional outline (or profile) of 
a Joukowski wing, we proceed as follows (Fig. 6). 

We first plot an xy system of coordinates with the ori- 
gin 0 such that the x axis forms the angle ^ with the hor- 
izontal direction of the wing and mark on the x axis the 
point L, for vmich x = - I, and on the y axis the point 
F, for which y = f . 

We now describe two circles and label them and K^. 
The center of the first circle is situated on the straight 

line LF at a distance 2 6 from the point F (beyond the 
section LF) . The circle, moreover, passes through the point 
L. The second circle likewise passes through the point L 
and its center is likewise on LF, the position of M 

on LF being determined by the following condition. If 0\ 
is the portion of the positive x axis cut off by the circle 
K^^ and- OVg the portion cut off by the circle K^, then 

OVi X OY2 f • 

We now draw, from ..the point 0, the two lines OA^ and 
OA2 , so as to form equal angles with' the x axis, A^^ being 
the point of intersection of the first lino with the circle 

and A s >the intersection of the second line vjit h- the circle 

* From "Zeitschrift fur Flugt echnik und Motorluf tschiff ahrt/* 
May 31, 1S13, pp. 130 and 131. 
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Kg. Then the center P of the line A-^Aq is the point sought 
on the- Joukowski v;ing profile- 
In plotting the preceding figures, 24 points v/ere found 
in this manner for each one, by shifting the first line from 
the point L 15*^ each time and drawing the second line sym- 
metrically with reference to the x axis. 

In order to determine the pressure on each point of the 
profile, when the wing is exposed to a horizontal wind having 
the velocity V, we must know the velocity q- at which the 
air flows oy each point of the profile. The pressure on each 
unit area of the wing surface is then proportional to • 

We can now find the values of q in a very simple man- 
ner- For this purpose, we draw a horizontal line through the 
point L. If we designate by h the distance of the point 
Ai (of the circle K^) from this horizontal lino, we obtain, 
for a^ny desired .point P of the figure, the corresponding 
value of q in the following manner. V/e take from the dia- 
gram the distance between the points A^ and A^ , at the mid- 
dle of which we had found the point P, and also the distances 
of the point A^ from the origin 0, from the center of 
the circle K^^ and from the horizontal line passing through 
L. We then have 

0 A, 3 h 
Ag >Ii Ai 

The mathematical proof for the given constructions is sim- 
ple. As already mentioned, the profile of a Joukowski v^ing 
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can be constructed by describing on the z .plane, with the aid 
of the formula z = ^ + -^-^ , the circle K, determined, by 
the camber and radii difference. Thia circle passes through 
the point ^ = - 

The systems of coordinates are plotted both in the t, 
plane and in the z plane in such manner that .the i axis and 
the X axis form the angle 3 with the horizontal wind direc- 
tion. 

If we now describe, in the s plane,, both circles, which 
we obtain from the given circle K in the plane by employ- 
ing the two conversion formulas 

f 

Zi = 2 ^ a.nd = ~ 

then these are the same two circles we designated above by 
Ki and Kg- 

The point Ai has the coordinate z^ and the point A3 
has the coordinate , hence the center of A;j A2 has the 

r 

coordinate z=-i(zi + Z2) = ^ + as desired. P is 

therefore an actual point on the Joukowski curve. 

The following formula holds good for the velocity q at 
which the air flows by every point on the Joukov/ski figure- 







) 




dz 











From z = t, + ly it follows that 
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whence v/e obtain 



A3A2 



0 Ai 

since the absolute value of - Zs equals the distance 

Ai A2 and the absolute value of .= the distance OAi . 

For ^ we obtain, from formula 2 of the preceding 

article, ^ = ^^ - ^ , in which h is the distance of the point 
M^Ai 

A^ from the horizontal line passing through L. In the ex- 
pression there given for the numerator, it is fequal to h and 
the denominator is equal to -^(M^Ai), as may be easily veri- 
fied. We thus obtain 

^ .M,A, A, A3 
which is just the formula given above for q. 



Translation by Dwight M. Miner, 
National Advisory Committee 
for Aeronautics* 
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Fig.l 



N.A.CA. Technical Memorandum No.336 Fig.3 




Fig. 2: f = 0, 6/i=L 
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